Decay of Dirac Massive Hair in the Background of Spherical Black Hole 
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The intermediate and late-time behaviour of massive Dirac hair in the static spherically general 
black hole spacetime was studied. It was revealed that the intermediate asymptotic pattern of decay 
of massive Dirac spinor hair is dependent on the mass of the field under consideration as well as the 
multiple number of the wave mode. The long-lived oscillatory tail observed at timelike infinity in 
the considered background decays slowly as i~ 5//6 . 

PACS numbers: 

I. INTRODUCTION 

Late-time behaviour of various fields in the spacetime of a collapsing body plays an important role in black hole 
physics, ft happens that regardless of details of the collapse or the structure and properties of the collapsing body 
the resultant black hole can be described only by few parameters such as mass, charge and angular momentum, black 
—i , holes have no hair. On its own, it is interesting to investigate how these hair loss proceeds dynamically. 
^ ' The research in this direction have quite long history. Price in [l[ showed that the late-time behavior is dominated by 
the factor i-( 2i + 3 ) : for each multipole moment I. It was revealed in Ref. Q that the decay-rate along null infinity and 
' along the future event horizon was governed by the power laws u~( l+2 ^ and v~ < -- l+3 \ where u and v were the outgoing 
Eddington-Finkelstein (ED) and ingoing ED coordinates. On the other hand, Ref. @ was devoted to the scalar 
perturbations on Reissner-Nordtrom (RN) background for the case when | Q \< M has the following dependence on 
time £-( 2i + 2 ) ; while for | Q \— M the late-time behavior at fixed r is governed by t~( l+2 \ 

A charged hair turned out to decay slower than a neutral one |4|-[6||. The late-time tails in gravitational collapse of a 
self-interacting (SI) fields in the background of Schwarzschild solution was reported by Burko 0] and in RN solution at 
intermediate late-time was considered in Ref. @ . At intermediate late-time for small mass m the decay was dominated 
by the oscillatory inverse power tails i~('+ 3 / 2 ) sin(mi). This analytic prediction was verified at intermediate times, 
where mM <C mt <C \/{mM) 2 . It was proved analytically [9| that for a nearly extreme RN spacetime the inverse 
power law behavior of the dominant asymptotic tail is of the form i~ 5 / 6 sin(mi). The asymptotic tail behaviour of 
SI scalar field was also studied in Schwarzschild spacetime [l(| ■ The oscillatory tail of scalar field has the decay rate 
of <~ 5 / 6 at asymptotically late time. The power-law tails in the evolution of a charged massless scalar field around 
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a fixed background of dilaton black hole was studied in Ref. [Ill], while the case of a self- interacting scalar field was 
elaborated in [l2| . The analytical proof of the above mentioned behaviour of massive scalar hair in the background 
of dilaton black hole in the theory with arbitrary coupling constant between U(l) gauge field and dilaton field was 
given in Ref. [22| . 

On its own, the late-time behaviour of massive Dirac fields were studied in the spacetime of Schwarzschild black hole 
[lH . It was found that the asymptotic behaviour of these fields is dominated by a decaying tail without any oscillation. 
The dumping exponent was independent on the multiple number of the wave mode and on the mass of the Dirac 
field. It also happened that the decay of the massive Dirac fields was slower comparing the decay of massive scalar 
field. The above analysis was supplemented by the studies of charged massive Dirac fields in the spacetime of RN 
black hole [l4| • The case of the decay of the fields in the stationary axisymmetric black hole background was studied 
numerically in Ref. [l5[ and it was found that in the case of Kerr black hole that the oscillatory inverse-power law 
of the dominant asymptotic tail behaviour is approximately depicted by the relation i -5 / 6 sin(mf). In Ref. |lq | both 
the intermediate late-time tail and the asymptotic behaviour of the charged massive Dirac fields in the background 
of Kerr-Newmann black hole was investigated. It was demonstrated that the intermediate late-time behaviour of the 
fields under consideration is dominated by an inverse power-law decaying tail without any oscillation. 



2 



The late-time behaviour of massive vector field obeying the Proca equation of motion in the background of 
Schwarzschild black hole was studied in jl?}. It was revealed that at intermediate late times, three functions char- 
acterizing the field have different decay law depending on the multiple number I. On the contrary, the late-time 
behaviour is independent on I, i.e., the late-time decay law is proportional to i -5 / 6 sin(mi). In Ref. [l8| the analytical 
studies concerning the intermediate and late-time decay pattern of massive Dirac hair on the dilaton black hole were 
conducted. Dilaton black hole constitutes a static spherically symmetric solution of the theory being the low-energy 
limit of the string theory with arbitrary coupling constant a. 

For the case of n-dimensional static black holes, the no-hair theorem is quite well established [l9j]. The mechanism 
of decaying black hole hair in higher dimensional static black hole case concerning the evolution of massless scalar 
field in the n-dimensional Schwarzshild spacetime was determined in Ref. (2p| . The late-time tails of massive scalar 
fields in the spacetime of n-dimensional static charged black hole was elaborated in Ref. [U] and it was found that 
the intermediate asymptotic behaviour of massive scalar field had the form £-( i +™/ 2 - 1 / 2 ). This pattern of decay was 
checked numerically for n = 5, 6. Quasi-normal modes for massless Dirac field in the background of n-dimensional 
Schwarzschild black hole were obtained in [23| , while in Ref. [HI studies of massless fermion excitations on a tensional 
threc-brane embedded in six-dimensional spacetime were provided. On the other hand, the intermediate and late-time 
decay of massive scalar hair on static brane black holes was elaborated in [25[ , while the decay of massive Dirac hair 
in the spacetime of black hole in question was studied in [26| . 

At the beginning of our universe several phase transitions leading to the topological defect formations might happened 
(27| . The interactions of topological defects such as global monopoles or cosmic strings with compact objects like 
black holes attract much attention and their physical characteristics are widely studied in literature. For instance, the 
black hole global monopole system has an unusual topological property of possessing a solid deficit angle [H|. The 
decay of massive scalar hair in the background of a Schwarzschild black hole with global monopole was considered in 
Ref. [29[ . It happened that the topological defects makes the massive scalar field hair decay faster in the intermediate 
regime comparing to the decay of such hair on the black holes without defects. On the other hand, the late-time 
behaviour was unaffected by the presence of global monopole. The Schwarzschild black hole global monopole system 
coupled to scalar fields was elaborated in Ref. (30j . 

Cosmic strings also acquire much interests as well as the cosmic string black hole systems. The metric of such a system 
assuming the distributional mass source was derived in [3l| (the so-called thin string limit). Then it was revealed 
that it constituted the limit of much more realistic situation when black hole was pierced by a Nielsen-Olesen vortex 
[32| . Especially interesting behaviour was found in the case of extremal black hole. Namely, for some range of black 
hole parameters it was shown that extremal black hole expelled the vortex (the so called Meissner effect). It turned 
our that in the case of extremal black holes in dilaton gravity one has always expulsion of the Higgs fields from their 
interiors [33| . The intermediate and late-time decay patterns of massive Dirac hair in the spacetime of black hole and 
topological defects were elaborated in Ref. [34| . The system of black hole with a global monopole and dilaton black 
hole pierced by a cosmic string was considered. Among all it was revealed that the late-time asymptotic decay of 
massive Dirac hair in the systems in question were proportional to i~ 5 / 6 and topological defects did not affect this 
relation. 

In our paper we shall discuss the intermediate and late-time behaviour of the massive Dirac spinor field in the 
spacetime of static spherically symmetric general kind of black hole. We revealed that the intermediate asymptotic 
behaviour is not the final pattern of the decay of the massive spinor hair. In Sec. II we present and summary our 
analytic arguments concerning the decay of massive spinor Dirac hair in the background of the black hole under 
consideration. In Sec. Ill we conclude that the long-lived oscillatory tail is generally observed at timelike infinity in 
static spherically symmetric black hole spacetime which decay pattern is proportional to . 

II. MASSIVE DIRAC HAIR ON SPHERICALLY SYMMETRIC BLACK HOLE 
A. Properties of the Dirac Equation 

The mass Dirac equation in a background of spherically symmetric black hole is given by 

^V^ip-mji(j = 0, (1) 

where is the covariant derivative V p = + juj^jajb, u and a are tangent and spacetime indices. There are 
related by e^, a basis of orthonormal one- forms. The quantity = oj ab are the associated connection one- forms 
satisfying de a + ojf, a A e b — 0, while the gamma matrices fulfill the relation {j a , j b } = 2n ab . Now we shall recall some 
basic properties of the Dirac operator^) = 7 M V M on an n-dimensional manifold (see e.g., [Uj]). In what follows we 



assume that the metric of the underlying spacetime may be rewritten as a product of the form 

g^dafda? = g ab (x)dx a dx b + g mn {y)dy m dy n . (2) 
The above metric decomposition will be subject to the direct sum of the Dirac operator, namely one obtains 

P=P x +P y . (3) 

If we define a Weyl conformal rcscaling as 

g M v = n 2 g^, (4) 

it consequently provides that one gets 

= n-i^+vpip , ^ = n-i^-^ijj. (5) 

Having in mind that spherically symmetric form of the line element provides also conformal flatness for a static metric, 
one obtains 

ds 2 = -A 2 dt 2 + B 2 dr 2 + C 2 dY? n ^ 2 , (6) 

where A — A{r), B — B(r), C = C(r) are functions only of the radial variable r, and the transverse metric dYi 2 l _ 2 is 
independent on t and on r. 

Suppose then, that ^ is a spinor eigenfunction on the (n — 2)-dimensional transverse manifold E. It leads to the 
following: 

p^ = A*. (7) 
In case of (n — 2)-dimensional sphere the eigenvalues for spinor 'J, where found in Ref . [35| . They imply 

A>=(< + ^) 2 , (8) 

where 1 = 0,1,... Using the properties given above one may also assume that the following is satisfied: 

pij) = mip. (9) 

It enables us to set the form of the spinor tp, i.e., 

i> = -JT—^TX®^- (10) 
A-2 c-r 1 

Next, the explicit calculations reveal that we arrive at the relation 

( 1 °d t + 1 1 d y )x = A(m- (11) 

where we have introduced the radial optical distance (i.e., the Regge- Wheeler radial coordinate) dy = B/Adr and 
7°,7 1 satisfy the Clifford algebra in two spacetime dimensions. 

We remark that an identical result may be obtained if a Yang-Mills gauge field A^ is present on the transverse manifold 
E. The only difference is that 



* = A^, (12) 



where p s A is the Dirac operator twisted by the the connection A^, Assuming that ijj cx e lut one obtains the second 
order equation of the form 

0W X = AV-^) 2 X. (13) 
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B. The Background of Spherically Symmetric Static Black Hole 



As was shown in Refs. [8j, [36j the spectral decomposition method will be fruitful for the analysis of the evolution of 
massive Dirac hair in the spacetime of static spherically symmetric black hole. The asymptotic tail is connected with 
the existence of a branch cut situated along the interval — m < u> < m. An oscillatory inverse power-law behaviour of 
the massive spinor field arises from the integral of Green function G(y, y'\ lu) around branch cut. The time evolution 
of the field may be written in the following form: 



X(V, t)= dy' 



G(y,y';t)My\0) + G t (y,y , - 1 t)4,(y',0) 



for t > 0, where the Green's function G(y, y'; t) is given by the relation 



dt 2 



dy 2 



V 



G(y,y'-t) = 5(t)5(y-y'). 



(14) 



(15) 



Next, our main task will be to find the black hole Green function. In the first step we reduce equation (|15p to an 
ordinary differential equation. To do it one can use the Fourier transform [3(| G(y, y'\ uf) = dt G(y, y'\ t)e lut . This 
Fourier's transform is well defined for Im u> > 0, while the corresponding inverse transform yields 



G(y,y';t) = — 



2tt 



oo+ze 



G(y,y';uj)e~ 



(16) 



oo+ie 



for some positive number e. The Fourier's component of the Green's function G{y,y';uj) can be written in terms of 
two linearly independent solutions for homogeneous equation as 



dy 2 



tS-V)xi = 0, i = 1,2, 



(17) 



where V = A 2 f m 



The boundary conditions for ipi are described by purely ingoing waves crossing the outer horizon H + of the n- 



dimensional static charged black hole ipi 

■02 ^ e -Vm 2 -^ 2 y at y . 



-iu>y 



as y 



-oo while 4>2 should be damped exponentially at i + , namely 



CO. 



In order to proceed further, we change the variables 



Xi 



A l/2 B -l/2 ' 



(18) 



where i = 1, 2. 

Although the formalism presented above is true for an arbitrary n-dimensional spacetime we shall consider general 
form of four-dimensional spherically symmetric static black hole. Black hole under consideration will be characterized 
by gravitational mass M and some other parameters M', Q, Q' of the background field in addition to the gravitational 
mass. 

We shall consider wave modes only in a far region from the black hole, as a generic behaviours observed in any black 
hole background. For that purpose, let us assume that r/M ^> 1, where M is the gravitational mass of background 
field. Having all these in mind, the expansion of the metric coefficients A(r) and B(r) as a power series in M/r give 
us the following: 



A M = 1 - — + ^2 



(19) 



while for the next metric function one has the relation 

2M' 



B{r) = 1 



(20) 



where by M' and Q, Q' we have denoted other parameters characterizing the background fields. 
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Let us assume that the observer and the initial data are situated far away from the considered black hole. We 
expand Eq. (|17[) as power series of M/r neglecting terms of order 0{uj/r) 2 ) and higher. It leads us to the following 
expression: 



where 



2a(ui 2 — to 2 ) + ol\ 



J 1 {a 2 + p 2 ) + Pi - m 2 a 2 - 2aa x - 2/3m 



ai = 2m(\ + Mm), Pi = A 2 + 4MtoA + Q 2 m 2 , 
a = M + M>, p = 91^1 + M> 2 -M 2 + MM , 



(21) 



(22) 



It turned out that the above equation can be solved in terms of Whittaker's functions, namely the two basic solutions 
are needed to construct the Green function, with the condition that | uj |> to. Consequently, it implies the result as 
follows: 



where we have denoted 



Xi = M K<jl (2Qr), X2 = W Ktjl (2ujr) 



jl = \Jl/A + u 2 (a 2 + P 2 ) + Pi-m 2 a 2 - 2aax - 2p~m 2 , 



(23) 



(24) 



2a(uj 2 — to 2 ) + ai 
2c5 



~ 2 2 2 

u) — m — uj . 



By virtue of the above relations the spectral Green function has the form 



G c (x,y;t) 



1 

2^ 



dw 



i r m 

— / dwf(Q) e-™\ 

^ J-rn 



W(w) 



(25) 



where W{Cj) is the Wronskian. 

First we discuss the intermediate asymptotic behaviour of the massive scalar field, i.e., when the range of parameters 
are M < r < i << M/{mM) 2 . The intermediate asymptotic contribution to the Green function integral gives the 
frequency equal to uj — 0{\Jm/t). It implies that k«1. One should have in mind that n stems from the 1/r term 
in the massive scalar field equation of motion. Thus it depicts the effect of backscattering off the spacetime curvature 
and in the case under consideration the backscattering is negligible. In the case of intermediate asymptotic behaviour 
one finally gets 



2 2 A- 1 r(-2 / 2) n\ + p) 
#r(2/i) r(| - p) 



(26) 



where we have used the fact that ujr <C 1 and the form of f(Q) can be approximated by means of the fact that 
M(a, 6, z) = 1 as z tends to zero. 

However, the resulting Green function can not be calculated analytically because of the fact that the parameter jx 
depends on uj. We calculate it numerically and results are provided in FigsQJS] In Fig. 1 we present the graph of the 
logarithm of intermediate Green function In | G c (r, r';t) \ versus time t for different mass parameter M = 0.1, M = 1.0 
and M — 5.0. In Fig. 2 we depict the dependence of the function in question on M' — 0.1, M' ~ 5.0, while in Fig. 3 
we get the relation among intermediate Green function and Q (for Q = 0.0 and Q = 1.0). Fig. 4 is connected 
with the dependence of the intermediate Green function on mass of the Dirac hair. We plotted the graph for 
to = 0.005, to = 0.01, to = 0.02. We finished with the Fig. 5, where the dependence on the multiple number of the 
wave mode was revealed for A = 0.01, A = 0.1 and A = 1.0. Summing it all up one can remark that the intermediate 
asymptotic Green function depends on the parameters M, M' , Q as well as the mass of the Dirac spinor field and the 
multiple number of the wave mode. 
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The different pattern of decay is expected when k ^> 1, for the late-time behaviour, when the backscattering off the 
curvature is important. In this case we use the limit 



M K , A (2u)r) m T(l + 2pt) (2uir)^K-i 1 J 2A (\/8Kwr). 

Consequently, /((D) yields 

r(i + 2/i)r(i-2A) , 

[rr ) 2 



2 ft 



J 2 ^(v / 8k^t) J_ 2 ^(-v/8KtDT-') — I 2 fi(V8Ktir) I-2fi(VSnu}r') 



(27) 



(28) 



(r(i + 2/2)) 2 r(-2/x) r(i + A-«) 

2/ir(2/i) r(i-/i-«;) 



(rr') 2 k 



-2A 



It can be noticed that the first part of the above Eq. ([28|) the late time tail is proportional to t 1 . Just we calculate 
the second term of the right-hand side of Eq. Q28p . For the case when n ^S> 1 it can be brought to the form written as 



G c ( 2 )(r,r';t) 



M 
2^ 



where we have defined 



while M provides the relation as follows: 

(r(i + 2/})) 2 r(-2£) 



M 



2/i T{2p) 



(rr'Y- 



1 + (_l)2A e -2«K 



J 2A (V8 



(29) 



(30) 



(31) 



At very late time both terms e llut and e 27rfc are rapidly oscillating. It means that the scalar waves are mixed states 
consisting of the states with multipole phases backscattered by spacetime curvature. Most of them cancel with each 
others which have the inverse phase. In such a case, one can find the value of G c (2) by means of the saddle point 
method. It could be found that the value 2ttk — wt is stationary at the value of w equal to the following: 



a 



it (2a(uj 2 - m 2 ) + ai) 



2V2rn 



Then approximating integration by the contribution from the very close nearby of ao is given by 



F 



im 2 / 3 V2 



(tt ai)i(mt)-i e lvia °\ 



(32) 



(33) 



In comparison to the late-time behaviour of the second term in Eq. (|28[) , the first term can be neglected. The dominant 
role plays the behaviour of the second term, i.e., the late-time behaviour is proportional to — |. Thus, the asymptotic 
late-time behaviour of the Green's function can be written in the form 



2\/2 m 2 / 3 5 i 5 . 

G c (r, r;t) = -= tt 6 a( (mt) 6 sm(mt) %(r, m) x(r , m). 

V3 



(34) 



The above relation provides the main conclusion of our investigations, i.e., analytical proof that the late-time asymp- 
totic decay pattern of massive Dirac field is unaffected by black hole characteristics and has the form of power law 
pattern proportional to £ -5 / 6 . 



III. CONCLUSIONS 



Dirac fields were intensively investigated in various contexts in spacetimes of black holes. Among all, problems con- 
cerning quasi-normal modes, high overtones were studied (see e.g., the latest works concerning this widely elaborated 



7 



field of research [38|]). In our paper we restrict our attention to the problem of the intermediate and late-time decay 
pattern of massive Dirac hair in the background of static spherically black hole. This problem is strictly bounded 
with the no-hair theorem for the objects in question. 

It has been shown in previous papers [i~3L llij that the oscillatory power-law decay rate of massive Dirac hair 
which decay rate is £~ 5 / 6 dominates at asymptotically late-time in Schwarzschild and RN black hole spacetimes. 
In our considerations we envisaged the fact that this behaviour was generic one observed in any static spherically 
symmetric background. We assume that a massive Dirac field propagates in a static spherically symmetric spacetime 
with asymptotically flat metric. The metric under consideration is described in terms of the ADM mass M and 
some other par ameters M', Q, Q' of the background fields in addition to the afore mentioned gravitational mass. 
In Refs.[l3l. Il4l. the complicated equation of motion for massive Dirac fields were solved by means of the so- 
called Newman-Penrose formalism. In our research we introduce much easier method of solving Dirac equation. The 
properties of the Dirac operator enable us to find the second order differential equation which in turn can be applied 
to construct the spectral Green function, crucial in our investigations. One should remark that our method may be 
used in arbitrary spacetime dimension contrary to the method of solvig Dirac Eqs. presented in the afore mentioned 
papers. However, we restrict our research to four-dimensional case. 

Unfortunately, the intermediate asymptotic behavior of massive spinor hair can not be calculated analytically. 
Then, we performed numerical calculations and found the graphs of the logarithm of the intermediate Green function 
versus time for changes of the parameters like M, M', Q, mass of the Dirac field m and the multiple number of the 
wave function. We revealed the dependence of the Green function in question on the above parameters. 

However, it turned out this is not the final pattern of vanishing the massive spinor hair. The resonance backscattering 
off the spacetime curvature reveals at very late times. It was proved analytically that the massive spinor decay in the 
spacetime under consideration is proportional to . So we conclude that this long-lived oscillating tail is generally 
observed at timclikc infinity in the black hole spacetimes. The same conclusion concerning the pattern of decay was 
achieved in Ref. [39| . where the massive scalar hair decay was elaborated in static spherically symmetric black hole 
spacetime. 

Our work comprises all the previously studied cases of static spherically symmetric black hole spacetimes on which 
propagates massive Dirac field. Finally, one should stress that we obtain the analytical proof that the late-time decay 
pattern of massive Dirac spinor hair in the background of general static spherically symmetric four-dimensional black 
hole it of the form t~ e . It does not depend on gravitational mass and additional parameters characterizing background 
fields. 

It will be interesting to study the decay pattern of higher spin fields in the general static spherically background (the 
research in this direction was begun in Ref. (l7| ) and to consider n-dimensional black hole background. We hope to 
return to these problems elsewhere. 
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FIG. 1: Amplitude of the Green's function vs. time for different Q. Parameters: M = 1, M' = 1, Q' = 0, rr' = 1000, A 
m = 0.01 



= 0.01, 




FIG. 2: Amplitude of the Green's function vs. time for different A. Parameters: M = 1, M' = 1, Q = 0, Q' = 0, rr' = 1000, 
m = 0.01. 
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FIG. 3: Amplitude of the Green's function vs. time for different M. Parameters: M' = 1, Q = 0, Q' = 0, rr' = 1000, A = 0.01, 
m = 0.01. 
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FIG. 4: Amplitude of the Green's function vs. time for different m. Parameters: M = 1, M' = 1, Q = 0, <Q' = 0, rr' = 1000, 
A = 0.01. 
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FIG. 5: Amplitude of the Green's function vs. time for different M' . Parameters: M = 1, Q = 0, Q' = 0, rr' = 1000, A = 0.01, 
m = 0.01. 



